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Abstract 

There is general agreement that the MiSaTaQuWa equations should describe the motion of a 
"small body" in general relativity, taking into account the leading order self-force effects. However, 
previous derivations of these equations have made a number of ad hoc assumptions and/or contain 
a number of unsatisfactory features. For example, all previous derivations have invoked, without 
O^' proper justification, the step of "Lorenz gauge relaxation" , wherein the linearized Einstein equation 

is written down in the form appropriate to the Lorenz gauge, but the Lorenz gauge condition is 
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then not imposed — thereby making the resulting equations for the metric perturbation inequiva- 
lent to the linearized Einstein equations. (Such a "relaxation" of the linearized Einstein equations 
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I is essential in order to avoid the conclusion that "point particles" move on geodesies.) In this 



paper, we analyze the issue of "particle motion" in general relativity in a systematic and rigorous 
way by considering a one-parameter family of metrics, gabW, corresponding to having a body 
(or black hole) that is "scaled down" to zero size and mass in an appropriate manner. We prove 
that the limiting worldline of such a one-parameter family must be a geodesic of the background 
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' metric, gab{^ = 0). Gravitational self- force — as well as the force due to coupling of the spin of 

the body to curvature — then arises as a first-order perturbative correction in A to this worldline. 
No assumptions are made in our analysis apart from the smoothness and limit properties of the 
one-parameter family of metrics, 5afe(A). Our approach should provide a framework for system- 
atically calculating higher order corrections to gravitational self-force, including higher multipole 
effects, although we do not attempt to go beyond first order calculations here. The status of the 
MiSaTaQuWa equations is explained. 
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I. INTRODUCTION 

The physical content of general relativity is contained in Einstein's equation, which has 
a well-posed initial value formulation (see, e.g., In principle, therefore, to determine 
the motion of bodies in general relativity — such as binary neutron stars or black holes — one 
simply needs to provide appropriate initial data (satisfying the constraint equations) on 
a spacelike slice and then evolve this data via Einstein's equation. However, in practice, 
it is generally impossible to find exact solutions of physical interest describing the motion 
of bodies by analytic methods. Although it now is possible to find solutions numerically 
in many cases of interest, it is difficult and cumbersome to do so, and one may overlook 
subtle effects and/or remain unenlightened about some basic general features of the solu- 
tions. Therefore, it is of considerable interest to develop methods that yield approximate 
descriptions of motion in some cases of interest. 

In general, the motion of a body of finite size will depend on the details of its composition 
as well as the details of its internal states of motion. Therefore, one can expect to get a 
simple description of motion only in some kind of "point particle limit". However, Einstein's 
equation is nonlinear, and a straightforward analysis shows that it does not make any 
mathematical sense to consider solutions of Einstein's equation with a distributional stress- 
energy tensor supported on a worldline^. Physically, if one tried to shrink a body down to 
zero radius at fixed mass, collapse to a black hole would occur before the point particle limit 
could be reached. 

Distributional stress-energy tensors supported on a world-line do make mathematical 
sense in the context of the linearized Einstein equation. Therefore, one might begin a treat- 
ment of gravitational self-force by considering a metric perturbation, hah-, in a background 
metric, Qab-, sourced by the stress-energy tensor of a "point particle" of mass M, given in 
coordinates (t,a;*) by 



Nevertheless, action principles corresponding to general relativity with point particle sources are commonly 
written down (see, e.g., eqs.(12.1.6) and (12.4.1)-(12.4.2) of [sj). There are no solutions to the equations of 
motion resulting from such action principles. By contrast, distributional solutions of Einstein's equation 
with support on a timelike hypersurface (i.e., "shells") do make mathematical sense 0, 
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where is the unit tangent (i.e., 4- velocity) of the worldhne 7 defined by = 
and r is the proper time along 7. (Here 5'^^^(x* — z'^(t)) is the "coordinate delta func- 
tion", i.e., j 5^^\x^ — z^{t))(Px^ = 1. The right side also could be written covariantly as 
8ttM J^64{x, z{r))ua{T)ub{r)dT where 64 is the covariant 4-dimensional delta-function and 
T denotes the proper time along 7.) However, this approach presents two major difficulties. 

First, the linearized Bianchi identity implies that the point particle stress-energy must 
be conserved. However, as we shall see explicitly in section HVl below, this requires that the 
worldline 7 of the particle be a geodesic of the background spacetime. Therefore, there are 
no solutions to equation ([T]) for non-geodesic source curves, making it hopeless to use this 
equation to derive corrections to geodesic motion. This difficulty has been circumvented 
m 

mm 

and other references by modifying ([T]) as follows. Choose the Lorenz gauge 
condition, so that equation ([1]) takes the form 

S^^^ir'' - zUf)) (]t 

V^Ve/lafe - 2i?U%d = -l67fMUa{t)Ub{t) ^—^=^ — , (2) 

V'hab = 0. (3) 

where hat = hat — \hgab with h = hahg""^- Equation ([2]) by itself has solutions for any source 
curve 7; it is only when the Lorenz gauge condition ([3]) is adjoined that the equations are 
equivalent to the linearized Einstein equation and geodesic motion is enforced. Therefore, if 
one solves the Lorenz-gauge form ([2]) of the linearized Einstein equation while simply ignoring 
the Lorenz gauge condition^ that was used to derive ([2]), one allows for the possibility non- 
geodesic motion. Of course, this "gauge relaxation" of the linearized Einstein equation 
produces an equation inequivalent to the original. However, because deviations from geodesic 
motion are expected to be small, the Lorenz gauge violation should likewise be small, and it 
thus has been argued Q| that solutions to the two systems should agree to sufficient accuracy. 

The second difficulty is that the solutions to eq. ([2]) are singular on the worldine of the 
particle. Therefore, naive attempts to compute corrections to the motion due to hab — 
such as demanding that the particle move on a geodesic of gab + hab — are virtually certain 
to encounter severe mathematical difficulties, analogous to the difficulties encountered in 



^ In some references, the failure to satisfy eq. ([3]) truly is ignored in the sense that it is not even pointed out 
that one has modified the linearized Einstein equation, and no attempt is made to justify this modification. 
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treatments of the electromagnetic self-force problem. 

Despite these difficulties, there is a general consensus that in the approximation that spin 
and higher multipole moments may be neglected, the motion of a sufficiently small body 
(with no "incoming radiation") should be described by self consistently solving eq. ([2]) via 
the retarded solution together with 

u'V.u'^ = -^{g'^'' + u^u''){2V,hlf - V,/if/)|^(^)«%^ , (4) 



where 



h'fix) = M - \9a^G\ (x, z{r')) u^' uUr' , (5) 



with G~^^^,y the retarded Green's function for eq. ([2]), normalized with a factor of — IGvr, 
following jol. The symbol r"^ indicates that the range of the integral extends just short 
of the retarded time Tret, so that only the "tail" (i-e., interior of the light cone) portion of 
the Green's function is used (see, e.g., reference j8[ for details). Equations ([2]) and (jlj) are 
known as the MiSaTaQuWa equations, and have been derived by a variety of approaches. 
However, there are difficulties with all of these approaches. 

One approach ^1 that has been taken is to parallel the analysis of [ol, [l^ in the electro- 
magnetic case and use conservation of effective gravitational stress-energy to determine the 
motion. However, this use of distributional sources at second-order in perturbation theory 
results in infinities that must be "regularized". Although these regularization procedures 
are relatively natural-looking, the mathematical status of such a derivation is unclear. 

Another approach [6] is to postulate certain properties that equations of gravitational self- 
force should satisfy. This yields a mathematically clean derivation of the self-force corrected 
equations of motion. However, as the authors of |d] emphasized, the motion of bodies in 
general relativity is fully described by Einstein's equation together with the field equations 
of the matter sources, so no additional postulates should be needed to obtain an equation of 
motion, beyond the "small body" assumption and other such approximations. The analysis 
given by jol shows that equation (jlj) follows from certain plausible assumptions. However, 
their derivation is thus only a plausibility argument for equation (jlj). Similar remarks apply 
to a later derivation by Poisson that uses a Green's function decomposition developed by 



Detweiler and Whiting 
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A third approach, taken by Mino, Sasaki, and Tanaka [sl and later Poisson — building 



on previous work of Burke 



Ul, d'Eath [12], Kates [IS;], Thorne and Hartle Q, and others 



involves the use of matched astymptotic expansions. Here one assumes a metric form in the 
"near zone" — where the metric is assumed to be that of the body, with a small correction due 
to the background spacetime — and in the "far zone" — where the metric is assumed to be that 
of the background spacetime, with a small correction due to the body. One then assumes 
that there is an overlap region of the body where both metric forms apply, and matches the 
expressions. The equations of motion of the body then arise from the matching conditions. 
However, as we shall indicate below, in addition to the "Lorenz gauge relaxation", there 
are a number of assumptions and steps in these derivations that have not been adequately 
justified. 

A more rigorous approach to deriving gravitational selfc^foree is suggested by tire work 
of Geroch and J ang [l5l | and later Geroch and Ehlers [IQ] on geodesic motion of small 



bodies (see also (17|). In [15|], one considers a fixed spacetime background metric gab and 
considers a smooth one-parameter family of stress-energy smooth tensors Tafe(A) that satisfy 
the dominant energy condition and have support on a world tube. As the parameter goes 
to zero, the world tube shrinks to a timelike curve. It is then proven that this timelike 
eurve must be a geodesie. Th,s result was general.zed in M to allow». to also vary with 
A so that Einstein's equation is satisfied. Within the framework of [16j, it therefore should 
be possible to derive perturbative corrections to geodesic motion, including gravitational 
self-force. However, the conditions imposed in IGj in effect require the mass of the body 
to go to zero faster than A^. Consequently, in this approach, a self- force correction like (jl]) 
to the motion of the body would arise at higher order than finite size effects and possibly 
other effects that would depend on the composition of the body. Thus, while the work of 
[igI provides a rigorous derivation of geodesic motion of a "small body" to lowest order, it is 
not a promising approach to derive gravitatational self-force corrections to geodesic motion. 



In this paper, we shall take an approach similar in spirit to that of 16|], but we will 
consider a different smooth, one-parameter family of metrics Qabi^), wherein, in effect, we 
have a body (or black hole) present that scales to zero size in a self-similar manner, with 



both the size and the mass of the body being proportional to A. In the 



the body (or black hole) shrinks down to a worldline, 7. As in [15|] and 16|, we prove that 



imit as A — i> 0, 



7 must be a geodesic of the "background spacetime" gab{^ = 0) — although our method of 
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proving this differs significantly from [15] and [16]. To first order in A, the correction to the 
motion is described by a vector field, Z*, on 7, which gives the "infinitesimal displacement" 
to the new worldline. We will show that, for any such one parameter family QabW, in the 
Lorenz gauge Z'^{t) satisfies 

(PZ^ 1 ■ ■ / ■ 1 A 



dt2 2M 

Here M and Sij are, respectively, the mass and spin of the body. The terms in parentheses 
on the right side of this equation correspond exactly to the gravitational self-force term in 
eq.(jl]); the first term is the Papapetrou spin- force [18]; the second term is simply the usual 
right hand side of the geodesic deviation equation. Equation ([6]) is "universal" in the sense 
that it holds for any one-parameter family satisfying our assumptions, so it holds equally 
well for a (sufficiently small) blob of ordinary matter or a (sufficiently small) black hole. 

Our derivation of ^ is closely related to the matched asymptotic expansions analyses. 
However, our derivation is a rigorous, perturbative result. In addition, we eliminate a num- 
ber of ad hoc, unjustified, and/or unnecssary assumptions made in previous approaches, 
including assumptions about the form of the "body metric" and its perturbations, assump- 
tions about rate of change of these quantities with time, the imposition of certain gauge 
conditions, the imposition of boundary conditions at the body, and, most importantly, the 
step of Lorenz gauge relaxation. Furthermore, in our approach, the notion of a "point 
particle" is a concept that is derived rather than assumed. It also will be manifest in our 
approach that the results hold for all bodies (or black holes) and that the physical results do 
not depend on a choice of gauge (although 2'*(r) itself is a gauge dependent quantity, i.e., the 
description of the corrections to particle motion depend on how the spacetimes at different 
A are identified — see section |V] and Appendix A). In particular, because the Lorenz gauge 
plays no preferred role in our derivation (aside from being a cornputationally convenient 



choice), our gauge transformation law is not, as in previous work [19|, restricted to gauges 
continuously related to the Lorenz gauge. Our approach also holds out the possibility of 
being extended so as to systematically take higher order corrections into account. However, 
we shall not attempt to undertake such an extension in this paper. 

Although ([6]) holds rigorously as a first-order perturbative correction to geodesic motion, 
this equation would not be expected to provide a good global-in-time approximation to 
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motion, since the small local corrections to the motion given by ([6]) may accumulate with 
time, and eq. IQ would not be expected to provide a good description of the perturbed 
motion when becomes large. We will argue in this paper that the MiSaTaQuWa equations, 
eqs. ([2]) and (jlj), should provide a much better global-in-time approximation to motion than 
eq.([n]), and they therefore should be used for self-consistent calculations of the motion of a 
small body, such as for calculations of extreme- mass-ratio inspiraF. 

We note in passing that, in contrast to Einstein's equation. Maxwell's equations are lin- 
ear, and it makes perfectly good mathematical sense to consider distributional solutions 
to Maxwell's equations with point particle sources. However, if the charge-current sources 
are not specified in advance but rather are determined by solving the matter equations of 
motion — which are assumed to be such that the total stress-energy of the matter and electro- 
magnetic field is conserved — then the full, coupled system of Maxwell's equations together 
with the equations of motion of the sources becomes nonlinear in the the electromagnetic 
field. Point particle sources do not make mathematical sense in this context either. It is 
for this reason that — despite more than a century of work on this problem — there is no 
mathematically legitimate derivation of the electromagnetic self-force on a charged parti- 
cle. The methods of this paper can be used to rigorously derive a perturbative expres- 
sion for electromagnetic self-force by considering suitable one-parameter families of coupled 
electromagnetic- matter systems, and we shall carry out this analysis elsewhere {20!]. How- 
ever, we shall restrict consideration in this paper to the gravitational case. 

This paper is organized as follows. In section [Tll we motivate the kind of one-parameter 
family of metrics, 5'afe(A), that we will consider by examining the one-parameter family of 
Schwarzschild-deSitter spacetimes with black hole mass equal to A. One way of taking the 
limit as A ^ yields deSitter spacetime. We refer to this way of taking the limit as the 
"ordinary limit". But we show that if we take the limit as A ^ in another way and 
also rescale the metric by A~^, we obtain Schwarzschild spacetime. We refer to this second 
way of taking the limit as A ^ as the "scaled limit". The scaled limit can be taken at 
any time to on the worldline 7. The simultaneous existence of both types of limits defines 
the kind of one parameter family of metrics we seek, wherein a body (or black hole) is 

^ This viewpoint is in contrast to that of [81] , where it is argued that second-order perturbations are needed 
for self-consistent evolution. 
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shrinking down to a world-line 7 in an asymptotically self-similar manner. The precise, 
general conditions we impose on QabW are formulated in section [Tll Some basic properties 
of Qabi^) that follow directly from our assumptions are given in section UTTl In particular, 
we obtain general "far zone" and "near zone" expansions and we show that, at any Iq, the 
scaled limit always yields a stationary, asymptotically fiat spacetime at A = 0. In section 
HVt we prove that 7 must be a geodesic of the "background spacetime" (i.e., the spacetime 
at A = resulting from taking the ordinary limit). In other words, to zeroth order in A a 
small body or black hole always moves on a geodesic. We also show that, to first order in A, 
the metric perturbation associated with such a body or black hole corresponds to that of a 
structureless "point particle". In section |Vl we define the motion of the body (or black hole) 
to first order in A by finding a coordinate displacement that makes the mass dipole moment 
of the stationary, asymptotically fiat spacetime appearing in the scaled limit vanish. (This 
can be interpreted as a displacement to the "center of mass" of the body.) In section IVll we 
then derive eq.(l6]) as the first order in A correction to 7 in the Lorenz gauge. (An appendix 
provides the transformation to other gauges.) Finally, in section IVlII we explain the status 
of the MiSaTaQuWa equation (jll). Our spacetime conventions are those of Wald ll, and we 
work in units where G = c = 1. Lower case Latin indices early in the alphabet (a, 6, c, ...) 
will denote abstract spacetime indices; Greek indices will denote coordinate components 
of these tensors; Latin indices from mid-alphabet {i,j,k,...) will denote spatial coordinate 
components. 



II. MOTIVATING EXAMPLE AND ASSUMPTIONS 

As discussed in the introduction, we seek a one-parameter family of metrics QabW that 
describes a material body or black hole that "shrinks down to zero size" in an asymptotically 
self-similar manner. In order to motivate the general conditions on Qabi^) that we shall 
impose, we consider here an extremely simple example of the type of one-parameter family 
that we seek. This example will provide an illustration of the two types of limits that we 
shall use to characterize gatW- 

Our example is built from Schwarzschild-deSitter spacetime, 

ds' = -(l~^- CoA dt'+fl-'^- Cy) dr' + r'dn' . (7) 
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(This metric, of course, is a solution to the vacuum Einstein's equation with a cosmological 
constant A = |Co rather than a solution with A = 0, but the field equations will not play 
any role in the considerations of this section; we prefer to use this example because of its 
simplicity and familiarity.) The desired one-parameter family is 

2^^^_ A 2MoA ^ o\ ,^2 , A 2MoA ^\ ^ 2 , „2.n2 



ds^X) = -[1- =^ - Cor'j dt' + (^1 - =^ - Cor'j dr' + r'd^\ (8) 

where we consider only the portion of the spacetime with r > XRq for some Rq > 2Mq. 
For each A, this spacetime describes the exterior gravitational field of a spherical body or 
black hole of mass AMq in an asymptotically deSitter spacetime. As A ^ 0, the body/black 
hole shrinks to zero size and mass. For A = 0, the spacetime is deSitter spacetime, which 
extends smoothly to the worldline r = 0, corresponding to where the body/black hole was 
just before it "disappeared'' 



As explained clearly in j2l|, the limit of a one-parameter family of metrics Qabi^) depends 
on how the spacetime manifolds at different values of A are identified. This identification 
of spacetime manifolds at different A can be specified by choosing coordinates for each 
A and identifying points labeled by the same value of the coordinates x^. If we use the 
coordinates (t, r, 6, (p) in which the one-parameter family of metrics ([8]) was presented to do 
the identification, then it is obvious that the limit as A — of Qabi^) is the deSitter metric 

ds\\ = 0) = - (1 - Cor^) dt" + (1 - Cyy^ dr^ + r^dn\ (9) 

This corresponds to the view that the body/black hole shrinks to zero size and mass and 
"disappears" as A ^ 0. 

However, there is another way of taking the limit of Qabi^) as A — > 0; the existence of 
this second limit is one of the key ideas in this paper. Choose an arbitrary time to and, for 
A > 0, introduce new time and radial coordinates hj i = (t — to)/A and f = r/X. In the new 
coordinates, the one-parameter family of metrics becomes 

ds\X) = -X'(l-^- CoX'fA dP + X'(l-^- CoX'fA ' df' + X'f'dn', f > Ro 



to- 
(10) 

We now consider the limit as A by identifying the spacetimes with different A at the 



9 



same values of the barred coordinates. It is clear by inspection of eq. (fTOi) that the the 
limit of Qabi^) as A ^ at fixed x'^ exists, but is zero. In essence, the spacetime interval 
between any two events labeled by Xi and goes to zero as A ^ because these events 
are converging to the same point on 7. Thus, this limit of QabW exists but is not very 
interesting. However, an interesting limit can be taken by considering a new one-parameter 
family of metrics (jabW by^ 

9fiu = >^~'^9fiu, (11) 

so that 

df{X) = (^1 - ^ - CoX'f'^ dP+(^l-^- CoX'f'^ ' df' + f'dn', f > Ro, (12) 
By inspection, the limit of this family of metrics is simply. 



ds\=o = - ( 1 - ^] dP+(l- ^] ' df' + f'dQ', f>Ro, (13) 



which is just the Schwarzschild metric with mass Mq. 

The meaning of this second limit can be understood as follows. As already discussed 
above, the one-parameter family of metrics ([8]) describes the exterior gravitational field of 
a spherical body or black hole that shrinks to zero size and mass as A ^ 0. The second 
limit corresponds to how this family of spacetimes looks to the family of observers placed 
at the events labeled by fixed values of x^. In going from, say, the A = 1 to the A = 1/100 
spacetime, each observer will see that the body/black hole has shrunk in size and mass by a 
factor of 100 and each observer also will find himself "closer to the origin" by this same factor 
of 100. Suppose now that this family of observers also "changes units" by this same factor 
of 100, i.e., they use centimeters rather than meters to measure distances. Then, except 
for small effects due to the deSitter background, the family of observers for the A = 1/100 
spacetime will report the same results (expressed in centimeters) as the observers for the 
A = 1 spacetime had reported (in meters). In the limit as A — > 0, these observers simply see 
a Schwarzschild black hole of mass Mq, since the effects of the deSitter background on what 
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A scaling of space (but not time) has previously been considered by Futamase [2 



context. Scaled coordinates also appear in the work of D'Eath 1^ and Kates [13 1 



in the post-Newtonian 
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FIG. 1: A spacetime diagram illustrating the type of one-parameter family we wish to consider, 
as well as the two limits we define. As A ^ 0, the body shrinks and finally disappears, leaving 
behind a smooth background spacetime with a preferred world-line, 7, picked out. The solid lines 
illustrate this "ordinary limit" of A ^ at fixed r, which is taken along paths that terminate away 
from 7 (i.e., r > 0). By contrast, the "scaled limit" as A — > 0, shown in dashed lines, is taken along 
paths at fixed f that converge to a point on 7. 

these observers will report disappear entirely in this limit. 

We will refer to the first type of limit (i.e., the limit as A ^ of QabW taken at fixed x^) 
as the ordinary limit of Qabi^), and we will refer to the second limit (i.e., the limit as A — > 
of A~^(7af,(A) taken at fixed x^) as the scaled limit of gabW- The simultaneous existence 
of both types of the above limits contains a great deal of relevant information about the 
one-parameter family of spacetimes ([8]). In essence, the existence of the first type of limit is 
telling us that the body/black hole is shrinking down to a worldline 7, with its mass going to 
zero (at least) as rapidly as its radius. The existence of the second type of limit is telling us 
that the body/black hole is shrinking to zero size in an asymptotically self-similar manner: 
In particular, its mass is proportional to its size, its shape is not changing, and it is not 
undergoing any (high frequency) oscillations in time. Figure [T] illustrates the two limits we 
consider. 

We wish now to abstract from the above example the general conditions to impose upon 
one-parameter families that express in a simple and precise way the condition that we have 
an arbitrary body (or black hole) that is shrinking to zero size — in an asymptotically self- 
similar way — in an arbitrary background spacetime. Most of the remainder of this paper will 
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be devoted to obtaining "equations of motion" for these bodies that are accurate enough 
to include gravitational self-force effects. A first attempt at specifying the type of one- 
parameter families gabW that we seek is as follows: 

• (i) Existence of the "ordinary limit": Qabi^) is such that there exists coordinates x° 
such that QnuiK x°') is jointly smooth in (A, x"), at least for r > for some constant 
R, where r = ^ERF = 1,2,3). For all A and for r > ^A, gab{\) is a vacuum 
solution of Einstein's equation. Furthermore, g^v{\ = 0, x°) is smooth in x°', including 
at r = 0, and, for A = 0, the curve 7 defined by r = is timelike. 

• (ii) Existence of the "scaled limit": For each to, we define i = (t ^ to)/\, x' = x'/A. 
Then the metric gfip{X;to;x") = X~^gjip{X;to;x"') is jointly smooth in (A, to; a;") for 
f = r/X > R. 

Here we have used the notation g^p to denote the components of g^b in the x" coordi- 
nates, whereas denotes the components of gab in the x" coordinates. It should be noted 
that, since the barred coordinates differ only by scale (and shift of time origin) from the 
corresponding unbarred coordinates, we have^ 

QfiD = X'^gi^„ . (14) 

Consequently, we have 

gupiX; to; t, x') = g^y{X; to + Xt, Xx') (15) 

since there is a cancelation of the factors of A resulting from the definition of gab and the 
coordinate rescahngs. It also should be noted that there is a redundancy in our description 
of the one-parameter family of metrics when taking the scaled limit: We define a scaled 
limit for all values of to, but these arise from a single one-parameter family of metrics gab{X)- 
Indeed, it is not difficult to see that we have 

^/ip(A; to; t + s, x') = gfiuiX; to + As; t, x') . (16) 

^ Note that in this equation and in other equations occurring later in this paper, we relate components of 
tensors in the barred coordinates to the corresponding components of tensors in unbarred coordinates. 
Thus, a bar appears over the indices on the left side of this equation, but not over the indices appearing 
on the right side of this equation. 



12 



In fact, our requirements on gab of the existence of both an "ordinary hmit" and a "scaled 
hnait" are not quite strong enough to properly specify the one-parameter families we seek. 
To explain this and obtain the desired strengthened condition, it is convenient to define the 
new variables 

a = r, /3 = \/r, (17) 

where the range of /? is < /3 < 1/-R. Let / denote a component of gabW in the coordinates 
x". However, instead of viewing / as a function of (A,x"), we view / as a function of 
{a, (3,t,d,(f)), where 6 and (p are defined in terms of by the usual formula for spherical 
polar angles. In terms of these new variables, taking the "ordinary limit" corresponds to 
letting /3 — * at any fixed a > 0, whereas taking the "scaled limit" corresponds to letting 



a — i> at any fixed f3 > (see figure 2(b)). Now, our assumptions concerning the ordinary 
limit imply that, at fixed {t, 9, (p) and at fixed a > 0, / depends smoothly on (3, including 
at /? = 0. On the other hand, our assumptions concerning the scaled limit imply that at 
fixed (t, 9, 0) and at fixed /? > 0, / is smooth in a. Furthermore, the last condition in the 
ordinary limit implies that for /5 = and fixed {t,9,(j)), f is smooth in a, including at a = 0. 

Thus, at fixed (t, 9, 0), our previously stated assumptions imply that / is well defined at 
the "origin" {a,P) = (0,0), and is smooth in a along the a-axis (i.e., /? = 0). However, our 
previously stated assumptions do not say anything about the continuity or smoothness of / 
as {a, (3) (0, 0) from directions other than along the a-axis. Such limiting behavior would 
correspond to letting r — > as A — > but at a rate slower than A, i.e., such that r/A — > oo. 
To see the meaning and relevance of this limiting behavior, let us return to our original 
motivating example, eq. ([8]) and take / to be the time-time component of this metric^. In 
terms of our new variables (fTTl) . we have 

f = -{l-2Mol3-Coa^) , (18) 



which is jointly smooth in (a,/?) at (0,0). By contrast, suppose we had a one-parameter 
family of metrics ^ab(A) that satisfies our above assumptions about the ordinary and scaled 

^ Note that if we wished to consider other components of this metric, we would have to transform back 
from the coordinates (r, 9^ 0) to "Cartesian-hke" coordinates that are well behaved at the origin = 
when A = 
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X B 

(a)The two limits in terms of r and A. A (b)The two limits in terms of a and /?. In the 

constant-A spacetime is shown as a thick line, new variables, the two types of limit appear on a 

The shaded region corresponds to the interior of more equal footing, with the ordinary and scaled 
the (shrinking) body, about which wc make no background metrics placed on cither axis. The 
assumptions. The ordinary limit is represented body is "pushed out" to finite (3, so that 

by solid lines and the scaled limit is represented assumptions made in a neighborhood of 

by dashed lines. While the ordinary background a = (3 = make no reference to the body, 
metric is on the r-axis, the scaled background 

metric is contained in the discontinuous behavior 
of the metric family at r = A = 0. 

FIG. 2: The two limits 

limits, but fails to be jointly smooth in (a, (3) at (0, 0). For example, suppose the time-time 
component of such a one-parameter family varied as 

In terms of the original variables (A, r), this corresponds to behavior of the form 

If we take the "ordinary limit" (A — > at fixed r > 0) of /, we find that / smoothly goes to 
— 1. Similarly, if we take the "scaled limit" (A ^ at fixed f — r/X> 0), we also find that 
/ smoothly goes to —1. However, suppose we let A — > but let r go to zero as r = cy/X. 
Then / will approach a different limit, namely c^/(l + c^). In essence, QabiX) corresponds 
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to a one-parameter family in which there is a "bump of curvature" at r oc V^. Although 
this "bump of curvature" does not register when one takes the ordinary or scaled limits, it 
is present in the one-parameter family of spacetimes and represents unacceptable limiting 
behavior as A — of this one-parameter family. 

In order to eliminate this kind of non- uniform behavior in A and r, we now impose the 
following addtional condition: 

• (iii) Uniformity condition: Each component of Qabi^) in the coordinates is a jointly 
smooth function of the variables {a, (3) at (0,0) (at fixed {t,6,(j))), where a and (3 are 
defined by eq.(fT7D. 

Assumptions (i)-(iii) constitute all of the conditions that we shall impose on (7afe(A). No 
additional assumptions will be made in this paper. 

We note that the coordinate freedom allowed by our conditions are precisely all coordinate 
t r ansf ormations 

x^" ^ x'^'i^x") (21) 

such that x'^(A, x*^) is jointly smooth in (A, x'^) for all r > CA for some constant C, and such 
that the Jacobian matrix dx'^/dx^ is jointly smooth in (a,/3) at (0,0) at fixed {t,6,(j)). 

It should be emphasized that our assumtions place absolutely no restrictions on the one- 
parameter family of spacetimes for r < XR, i.e., this portion of these spacetimes could 
equally well be "filled in" with ordinary matter or a black hole^. It also should be noted 
that the "large r" region of the spacetime will not be relevant to any of our considerations, 
so it is only necessary that our conditions hold for r < K for some constant K. 

Finally, although it may not be obvious upon first reading, we note that our assump- 
tions concerning QabW are closely related to the assumptions made in matched asymptotic 
expansion analyses. As we shall see in the next section, in essence, our assumption about 
the existence of an ordinary limit of gab{X) corresponds to assuming the existence of a "far 

Indeed, it could also be "filled in" with "exotic matter" (failing to satisfy, say, the dominant energy 
condition) or a naked singularity (of positive or negative mass). However, in cases where there fails to be 
a well posed initial value formulation (as would occur with certain types of exotic matter and with naked 
singularities) and/or there exist instabilities (as would occur with other types of exotic matter), if is far 
from clear that one should expect there to exist a one-parameter family of solutions gabW satisfying our 
assumptions. 
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zone" expansion; our assumption about the existence of a scaled limit of gabW corresponds 
to assuming the existence of a "near zone" expansion; and our uniformity assumption corre- 
sponds closely to the assumption of the existence of a "buffer zone" where both expansions 
are valid. 

III. CONSEQUENCES OF OUR ASSUMPTIONS 

In this section, we derive some immediate consequences of the assumptions of the previous 
section that will play a key role in our analysis. These results will follow directly from the 
"uniformity condition" and the consistency relation (|T6|) . 

Since, by the uniformity assumption, the coordinate components of the one-parameter 
family of metrics Qabi^) are jointly smooth in the variables [a, (3) at (0, 0), we may approx- 
imate g^y by a joint Taylor series in a and (3 to any finite orders and M by 

N M 

g,y{\- 1, r,e,<p) = Y,Yl (^y^Mnmit, e, 0) + 0(a^+i) + 0(/3^+^) . (22) 

n=0 m=0 

Substituting for a and (3 from eq. fll7l) . we have 

(7^,(A;t,r,0,(/.) = ^^r"(-j (a^.)„^(t, 0, </>) , (23) 

n=0 m=0 ^ ^ 

where here and in the following, we drop the error term. We can rewrite this equation as a 
perturbation expansion in A: 

M N 

g,y{X;t,r,e,<P) = ^ A'" r"-'"(a^,)„„(t, ^, 0) . (24) 

m=0 n=0 

We will refer to eq. flM|) as the far zone expansion of gab{X)- It should be noted that the 
mth-order term in A in the far zone perturbation expansion has leading order behavior of 
l/r"^ at small r. However, arbitrarily high positive powers of r may occur at each order 
in A. Finally, we note that the angular dependence of {a^y)no(t, 9, 0) is further restricted 
by the requirement that the metric components g^u{X = 0) must be smooth at r = when 
re-expressed as functions of x\ In particular, this implies that (a^i/)oo cannot have any 
angular dependence. 
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Equivalently, in view of eq. (fT5!) . we can expand g^ip as 



N M 



ivi / 1 \ rn—n 

= J] J] A" ( - j {a^u)nm{to + Xt, e, 0) . (25) 

n=0 m=0 ^ ^ 

By further expanding {a^^)nm in ^ about t = 0, we obtain 

NMP . s m-n 

-g-,,{\- to; i, r,e,<P) = J2J2Yl ( ) ( V)nn^p(io, 0, 0) , (26) 

n^n m=n n=n \ / 



n=0 m=0 p=0 

where 



1 9P 



(27) 



t=to 

We can rewrite this as a perturbation series expansion in A: 



N+P min(g,P) M 

/ -I \ rn—q+p 

^^p(A;to;t,r,^,0)= 5^ A^ 5Z ^ ) (V)(9-pw(^o; ^, 0) • (28) 

5=0 p=0 m=0 ^ ^ 



We will refer to eg. (1281) as the near zone expansion of QabW- We see from this formula that 
the scaled metric, viewed as a perturbation series in A, follows the rule that the combined 
powers of i and f are allowed to be only as positive as the order in perturbation theory. By 
contrast inverse powers of f of arbitrarily high order are always allowed. Of course, only 
non-negative powers of i can occur. 

By inspection of eq. fl28l) . we see that the "background" (A = 0) scaled metric is given by 



M 

I I \ 

a^i/)om(to;6',0) , (29) 



m=0 

where we have used the fact that (6^i/)omo = (ct^ti/)om- Thus, we see that there is no de- 
pendence of gp.i,{\ = 0;to) on t and only non-positive powers of f occur. Thus, we see that 
(jabi^ = 0) is a stationary, asymptotically flat spacetime. However, the limiting, stationary, 
asymtotically flat spacetime that we obtain may depend on the choice of the time, to, on 
the worldline, 7, about which the scaling is done. 

Our "far zone expansion" , eg. (1211) . appears to correspond closely to the far zone expansion 
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used in matched asymptotic expansion analyses [5|, |8|]. However, our "near zone expansion" 
differs in that we define a separate expansion for each time to rather than attempting a 
uniform in time approximation with a single expansion. Such expansions require an ad- 
ditional "quasi-static" or slow-time variation assumption for the evolution of the metric 
perturbations. A further difference is that the conclusion that the background (A = 0) met- 
ric is stationary and asymptotically fiat has been derived here rather than assumed. Indeed, 
in other analyses, a particular form of the background metric (such as the Schwarzschild 
metric) is assumed, and the possibility that this metric form might change with time (i.e., 
depend upon to) is not considered. In addition, in other analyses boundary conditions at 
small f (such as regularity on the Schwarzschild horizon) are imposed. In our analysis, we 
make no assumptions other than the assumptions (i)-(iii) stated in the previous section. In 
particular, since we make no assumptions about the form of the metric for f < R, we do not 
impose any boundary conditions at small f. 

Finally, it is also useful to express the consistency relation (fT6l) in a simple, differential 
form. We define 

iKf,^)npmiX;to;t,x') = i—j {-g^j {g^j 9i>^; ^0] t, x')\x=0,to,i=0 (30) 

Then, a short calculation shows that 

K(n+l)p{m+l) = in+ l)Kn(p+l)m (31) 

as well as 

Knpm = if n < m . (32) 

Setting 77, = 0, we see that the last relation implies that g{X = 0; to; t, a;*) is stationary, as we 
have already noted. 



IV. GEODESIC MOTION 



In this section, we will prove that the worldline 7 appearing in assumption (i) of section 
ininiust, in fact, be a geodesic of the background metric gab{^ = 0). This can be interpreted 
as establishing that, to zeroth order in A, any body (or black hole) moves on a geodesic 
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of the background spacetime. In fact, we will show considerably more than this: We will 
show that, to first order in A, the far zone description of gab{^) is that of a "point particle". 
As previously mentioned in the Introduction, our derivation of geodesic motion is similar in 
spirit to that of [16] in that we consider one-parameter families of solutions with a worldtube 
that shrinks down to a curve 7, but the nature of the one-parameter families that we consider 



here are very different from those considered by 16|], and our proof of geodesic motion is 



very different as well. Our derivation of geodesic motion also appears to differ significantly 
from pervious derivations using matched asymptotic expansions j^, [s], Q, 3]. 

We begin by writing the lowest order terms in the far zone expansion, eq. (12^ . as follows: 



QafsW = (aa/3)oo(i) + (aa/3)io(i, 0, (j))r + 0(r 



2^ 



+ A 



(aa/3)oi(i, + {aap)ii{t, 9, (p) + 0(r) 



+ O(A2),r>0, (33) 



where we have used the fact that (00/3)00 can depend only on t, as noted in the previous 
section. Since the worldline 7, given by = 0, was assumed to be timelike* in the spacetime 
Qabi.^ = 0), without loss of generality, we may choose our coordinates so that Qapi,^ = 
0, X* = 0) = (aa/3)oo(^) = Vai3- (One such possible choice of coordinates would be Fermi 
normal coordinates with respect to 7 in the metric Qabi^ = 0). We emphasize that we make 
the coordinate choice {aa/3)oo = f]af3 purely for convenience — so that, e.g., coordinate time 
coincides with proper time on 7 — but it plays no essential role in our arguments.) Choosing 
these coordinates, and letting ha/3 denote the 0(A) piece of the metric, we see that ga/sW 
takes the form 

9ap = Vap + 0{r) + \Kp + 0{\^) , (34) 



where 



Cal3{t,e,(j)) 

ha/3 = h 0(1) . (35) 

r 



where in eq.f l35p . the term "0(1)" denotes a term that, when multiplied by for any e > 0, 
vanishes as r ^ 0. 

Now, by assumption (i) of section [Tll for each A, gabW is a vacuum solution of Einstein's 



We made this assumption explicitly in condition (i) of section |TT1 However, if, instead, we had assumed 
that the "interior region" r < XR, were "filled in" with matter satisfying the dominant energy condition, 
then it should be possible to prove that 7 must be timelike. 
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equation for r > XR and is jointly smooth in (A, x") in this coordinate range. It follows that 
for all r > 0, hat is a solution the linearized Einstein equation off of gab{X = 0), i.e., 

G^blhcd] = —VaVth', - ^Wchab + W^tha)c = 0, f > , (36) 

where here and in the following, Va denotes the derivative operator associated with gab{X = 
0), and indices are raised and lowered with gab{X = 0). Equation (!36|) holds only for r > 0, 
and, indeed, if Cai3 7^ 0, hab is singular at r = 0. However, even if Cais 7^ 0, the singularity of 
each component of hab is locally integrable with respect to the volume element associated 
with QabiX = 0), i.e., each component, hajs, is a locally function on the entire spacetime 
manifold, including r = 0. Thus, hab is well defined as a distribution on all of spacetime. The 
quantity Tab = Gl^blhcd]/^^ is therefore automatically well defined as a distribution. This 
quantity has the interpretation of being the "source" for the metric perturbation fl35l) — even 
though all of our spacetimes gab{X) for A > have excluded the "source region" r < A^. It 
follows immediately from eq.( l36i) that, as a distribution. Tab must have support on 7 in the 
sense that it must vanish when acting on any test tensor field whose support does not 
intersect 7. We now compute Tab- 
By definition. Tab = G^ab[hcd]/8Tc is the distribution on spacetime whose action on an 
arbitrary smooth, compact support, symmetric tensor field fed = fdc is given by 



^^nf) 



[ Gi;)[/,d/i"^v^rf^x = 0, (37) 



where ^J—gd^x denotes the volume element associated with gab{X = 0) and we have used 
the fact that the operator G^^ij is self- adjoint^. Note that the right side of this equation is 
well defined since G^^fj[fcd\ is a smooth tensor field of compact support and h"-^ is locally 
L^. We can evaluate the right side of eq. fl57|) by integrating over the region r > e > and 
then taking the limit as e — 0. In the region r > e, h"-^ is smooth, and a straightforward 
"integration by parts" calculation shows that 

Gll^l/cJ/^-'-G^Mr = V,.% (38) 



23 | for the definition of adjoint being used here. If G^^^ were not self-adjoint, then the adjoint of G^^^ 



See 

would have appeared in eq. ([57|) . 
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where 



+ V^/i/ - W7 + hVbf' - Vkhf' , (39) 
where / = fabg"'^{^ = 0). Since G^^^[hcd\ = for r > 0, it follows immediately that 

Using eqs. (I35|) and ( l39l) . we find that T(/) takes the form 

r(/)= J dtNab{t)r\t,r = 0) , (41) 

where Nab{t) is a smooth, symmetric {Nab = N^a) tensor field on 7 whose components are 
given in terms of suitable angular averages of Cajj and its first angular derivatives. In other 
words, the distribution Tat is given by^° 

Tab = Nabit)—^— , (42) 
■\/ g clt 

where 5^^\x'') is the "coordinate delta-function" (i.e., J 5^^\x'')d'^x^ = 1). 

We now use the fact that, since the differential operator G^^i^ satisfies the linearized 
Bianchi identity V^G^^ij = 0, the distribution Tab must satisfy V^Tab = in the distributional 
sense. This means that the action of Tab must vanish on any test tensor field of the form 
fab = V(a/b) where is smooth and of compact support. In other words, by eq. fKT]) . the 
tensor field Nab on 7 must be such that for an arbitrary smooth vector field /" on spacetime, 
we have 

J dtNabm''f\t,r = 0)=0 . (43) 

Now for any i = 1,2, 3, choose /" to have components of the form = x^F{x^, x^, x^)c'^{t), 
where each c'^ (/i = 0, 1, 2, 3) is an arbitrary smooth function of compact support in t and F 
is a smooth function of compact spatial support, with F = 1 in a neighborhood of 7. Then 



In fact, by our coordinate choice, we have ^—g = 1 on 7 and ^ — 1, but we prefer to leave in these 



dt 

factors so that this formula holds for an arbitrary choice of coordinates. 
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eq. (03]) yields 

j dtN,^{t)c^{t) = (44) 

for all c^(t), which immediately implies that A^j^ = N^i = for all i = 1,2,3 and all 
/i = 0, 1, 2, 3. In other words, we have shown that Nab{t) must take the form 

Nab = M{t)UaUb (45) 

where denotes the unit tangent to 7, i.e., u"" is the 4- velocity of 7. Now choose to be 
an arbitrary smooth vector field of compact support. Then eqs. fj43l) and fj45l) yield 

= J dtM{t)Ub{UaV''f) =- J dtu''Va{M{t)Ub)f , (46) 

where we integrated by parts in t to obtain the last equality. Since is arbitrary, this 
immediately implies that 

u''Wa{M{t)ub) = . (47) 

This, in turn, implies that 

dM/dt = , (48) 
i.e., M is a constant along 7, and, if M 7^ 0, 

u'^Vau' = , (49) 

i.e., in the case where M 7^ 0, 7 is a geodesic of Qabi^ = 0), as we desired to show 

In summary, we have shown that for any one-parameter family of metrics gabW satisfying 
assumptions (i)-(iii) of section [TTl to first order in A, the far zone metric perturbation hab 
corresponds to a solution to the linearized Einstein equation with a point particle source 

Tab = MuaUb^^— , 50 
^J—g dt 



Some previous derivations of geodesic motion do not appear to make explicit use of the fact 

that M 7^ 0. It is critical that this assumption be used in any valid derivation of geodesic motion, since a 
derivation that holds for Af = effectively would show that all curves are geodesies. 
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where M is a constant and u"" is the 4-velocity of 7, which must be a geodesic if M 7^ 0. 
We refer to M as the mass of the particle. It is rather remarkable that the point particle 
source (1501) is an output of our analysis rather than an input. Indeed, we maintain that 
the result we have just derived is what provides the justification for the notion of "point 
particles" — a notion that has played a dominant role in classical physics for more than three 
centuries. In fact, the notion of point particles makes no mathematical sense in the context 
of nonlinear field theories like general relativity. Nevertheless, we have just shown that the 
notion of a (structureless) "point particle" arises naturally as an approximate description of 
sufficiently small bodies — namely, a description that is valid to first order in A in the far zone 
for arbitrary one-parameter families of metrics QabW satisfying the assumptions of section 
[m This description is valid independently of the nature of the "body", e.g., it holds with 
equal validity for a small black hole as for a small blob of ordinary matter. 



V. DESCRIPTION OF MOTION TO FIRST ORDER IN A 



In the previous section, we established that, to zeroth order in A, any body (or black 
hole) of nonvanishing mass moves on a geodesic of the background spacetime. Much of the 
remainder of this paper will be devoted to finding the corrections to this motion, valid to 
first order in A in the far zone. In this section, we address the issue of what is meant by the 
"motion of the body" to first order in A. 

The first point that should be clearly recognized is that it is far from obvious how to 
describe "motion" in terms of a worldline for A > 0. Indeed, the metric QabW is defined only 
for r > XR, so at finite A the spacetime metric may not even be defined in a neighborhood 
of 7. If we were to assume that R» M and that the region r < A^ were "filled in" with 
sufficiently "weak field matter" — so that TZR^ « 1, where TZ denotes the supremum of 
the components of the Riemann curvature tensor of gab{X) in the "filled in" region — then it 
should be possible to define a "center of mass" worldline at finite A, and we could use this 
worldline to characterize the motion of the body 2J] . However, we do not wish to make any 
"weak field" assumptions here, since we wish to describe to motion of small black holes and 
other "strong field" objects. Since it is not clear how to associate a worldline to the body 
at finite A — and, indeed, the "body" is excluded from the spacetime region we consider at 
finite A — it is not clear what one would mean by a "perturbative correction" to 7 to first or 
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higher order in A. 

A second point that should be understood is that if we have succeeded in defining the 
worldhnes describing the motion the body at finite A, 

x\X, t) = z\X, t) = \Z\t) + 0(A2) , (51) 

then the "first order in A perturbative correction", Z*, to the zeroth order motion 7 (given 
by x*(t) = 0) is most properly viewed as the spatial components of a vector field, Z", defined 
along 7. This vector field describes the "infinitesimal displacement" to the corrected motion 
to first order in A. The time component, of Z"- depends on on how we identify the 
time parameter of the worldlines at different values of A and is not physically relevant; we 
will set Z'^ = so that Z'^ is orthogonal to the tangent, m'*, to 7 in the background metric 
Qabi^ = 0). Thus, when we seek equations of motion describing the first order perturbative 
corrections to geodesic motion, we are seeking equations satisfied by the vector field Z°'(t) 
on 7. 

A third point that should be clearly recognized is that Z'^ and any equations of motion 
satisfied by Z° will depend on our choice of gauge for hab- To see this explicitly, suppose 
that we perform a smooth^^ gauge transformation of the form 

^ = _ AA^(x'^) + 0(A2) . (52) 



Under this gauge transformation, we have 

hf,u = h^^ - 2V(^A^) . (53) 

However, clearly, the new description of motion will be of the form [l^ 

x\\,t)=z\\t) (54) 



Gauge transformations where is not smooth at a;* = are also permitted under the coordinate freedom 
stated at the end of sectionlTTl However, it suffices to consider smooth for our consderations here. The 
change in the description of motion under non-smooth gauge transformations will be treated in Appendix 
A. 
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where 

z\t) = z\t) - \A\t, = 0) + 0{\^) . (55) 
Thus, we see that Z"" transforms as 

Z\t) ^ Z\t) = Z\t) - A\t, = 0) (56) 

in order that it describe the same perturbed motion. Thus, the first order correction, Z°'{t), 
to the background geodesic motion contains no meaningful information by itself and, indeed, 
it can always be set to zero by a smooth gauge transformation. Only the pair [hat, Z"'{t)) 
has gauge invariant meaning. 

We turn now to the definition of the first order perturbed motion. Our definition relies 
on the fact, proven in section UTTl above, that for each Iq, gfip{X = 0;to;x°') is a stationary. 



asymptotically fiat spacetime. Therefore, gabi^ = 0) has well defined sets o: 
parity") and angular momentum ("magnetic parity") multipole moments 



mass ("electric 
261 ] and, in- 
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deec 



271, 



, the spacetime is characterized by the values of these two sets of multipole moments 



281]. The multipole moments (other than the lowest nonvanishing multipoles of each 
type) depend upon a choice of conformal factor [25|, l26|], which, rougly speaking,corresponds 
to a choice of "origin". We choose the conformal factor f2 = to define all of the multi- 
poles, corresponding to choosing the origin at f = 0. For a metric of the form eq. (1291) — with 
(o/i!/)oo = Vfif by our coordinate choice imposed in the previous section that g^^i^ = 0) = 1]^^, 
on 7 — the mass will be simply the / = part of the coefficient of 1/f in the large f expansion 
of |(?ft(A = 0; to)- Similarly, the mass dipole moment will be the coefficient of the / = 1 part 
of this quantity at order 

It is well known that if the mass is nonzero, the mass dipole moment is "pure gauge" 
and can be set to zero by choice of conformal factor/ "origin". We now explicitly show that, 
with our choice of conformal factor Q = the mass dipole moment can be set to zero 

by a smooth gauge transformation of the form (!52l) . It follows from the linearized Einstein 
equation fl36|) with source fl50|) applied to hab, eq. fl35|) . that the time-time component of hab 
takes the form 

hu='^ + 0(1) , (57) 
i.e., in the notation of eq. flM|) . we have {att)oi = 2M. Comparing with eq. fl28|) (and also 
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using the fact that {att)oo = —1), we see that at each to 

2M 



9h{X = 0; to) = -(1 - — ) + 0{l/n . (58) 

From this equation, we see that the "particle mass" , M, of the "source" of the far zone metric 
perturbation (see eq. fl50|) ) is also the Komar/ADM mass of the stationary, asymptotically 
flat spacetime (jabi^ = 0;to). We now calculate the effect of the coordinate transforma- 
tion ( !52l) on cjui^ = 0;to)- The transformation ( l52l) corresponds to changing the barred 
coordinates by 

xf' ^^^^ = x^'- A^ito, x' = 0) + 0(A) , (59) 
i.e., to zeroth order in A, it corresponds to a "constant displacement" of coordinates. Since 

111 AiX' 



r \x' + A'{to,0)\ f f3 



+ 0(1/P) , (60) 



it can be seen that this transformation have the effect of changing the mass dipole moment 
by — My4*. In particular, we can always choose so as to set the mass dipole moment to 
zero. 

Now, the "near zone" coordinates x^ for which the mass dipole moment vanishes have 
the interpretation of being "body centered" coordinates to zeroth order in A. The origin 

= of the corresponding "far zone" coordinates therefore has the interpretation of 
representing the "position" of the center of mass of the body to first order in A. We shall 
use this to define the correction to geodesic motion to first order in A by proceeding as 
follows: 

First, we shall choose our coordinates, x^, to zeroth order in A by choosing convenient 
coordinates for the "background spacetime" Qabi^ = 0). (We will use Fermi normal coordi- 
nates based on 7.) Next, we will define our coordinates, x^, to first order in A by choosing a 
convenient gauge for hat, eq. fl35l) . (We will choose the Lorenz gauge V"(/?.a{, — \hgab) = 0.) 
Then we will introduce the smooth coordinate transformation ( |52l) . and impose the require- 
ment that A^ be such that the mass dipole moment of gab{^ = 0; to) vanish for all to. Since 
the "location" of the body in the new coordinates is z*(t) = 0, the first order perturbative 
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correction Z°'{t) to the motion of the body in our original coordinates will be given by 

Z\t) = A\t,x^ = 0) (61) 

Of course, the particular Z°'{t) that we obtain in any given case will depend upon the 
particular one-parameter family Qabi^) that we consider. What is of interest is any "universal 
relations" satisfied by ^'"(t) that are independent of the choice of one-parameter family 
satisfying assumptions (i)-(iii) of section [Tll Such universal relations would provide us with 
"laws of motion" for point particles that take self-force effects into account. In the next 
section, we will show (via a lengthy calculation) that such a universal relation exists for 
(Pzydt"^, thus providing us with general "equations of motion" for all "point particles", 
valid to first order in A. 

Finally, we note that if we wish to describe motion beyond first order in A, it will be 
necessary to define a "representative world line" in the far zone to at least second order in 
A. We shall not attempt to do so in this paper. The definition of a suitable representative 
worldline is probably the greatest obstacle to extending the results of this paper to higher 
order in perturbation theory. 

VI. COMPUTATION OF PERTURBED MOTION 

In the section IIVI we found that first-order far zone perturbations of the background 
spacetime gab{^ = 0) are sourced by a point particle stress-energy, eq. (l50l) . For the re- 
mainder of this paper, we will assume that M 7^ 0, so that, as shown in section IIVI the 
lowest order motion is described by a geodesic, 7, of the background spacetime. We will 
need expressions for the components of the far zone metric, g^u\x=o, its first order pertur- 
bation, hf^i, = dg^u/d\\\=Q, and its second order perturbation j^^ = |c^^fi'^i//<9A^|A=o- It is 
convenient to choose our coordinates to zeroth order in A to be Fermi normal coordinates 
with respect to the background geodesic 7, and choose these coordinates to first order in A 
so that hfj_y satisfies the Lorenz gauge condition V^hab = 0, where hab = hab — ^hgab\\=o with 
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h = h 

abg"'''\x=o- Then the hnearized Einstein equation reads 



V''VXb-2R\Ad = -l6nM I 5i{x,z{r))ua{T)ui,{T)dT, (62) 
V'hab = 0. (63) 



This system of equations can be solved using the Hadamard expansion techniques of 
DeWitt and Brehme p, [g], 1^. Since this technology has been used in all previous derivations 
of gravitational self-force, we do not review it here but simply present results. Equation 
(2.27) of Mino, Sasaki, and Tanaka [5] provides a covariant expression for the perturbations in 
terms of parallel propagators and Synge's world function on the background metric (see, e.g., 
reference |8i] for definitions of these quantities). The Fermi normal coordinate components 
of these tensors are easily calculated with the aid of expressions from section 8 of Poisson 
js]. Combining this with the form of the background metric in Fermi normal coordinates, 
we obtain 

5fa/3(A; t, X") = + Baipj{t)x'- x^ + O(r^) 



+ A + h^^it, 0) + h%{t, 0)x* + MTZ^pit, x') + O(r^) J + O(A^) , (64) 

where the quantities B^p-ys and TZap are defined by the following expressions in terms of the 
Fermi normal coordinate components of the Riemann tensor of the background metric 

x^x^ 

Bom = —Rom "^oo = 7 Rom— (65) 
Bim = —-^^Rim TZio = —Rim "^Riokox^ (66) 



3 
1 





X^X^ 


"^00 


— 7Rom 




r 




2 x^x^ 


TZio 


— ^Rim 




3 r 




- ^ 


2 ^ikjl 



Bikji — —-Rikji '^ij ~ ~~-^ikji~ 4:rRiQjo, (67) 



and /i^^ and /i^^^ are given by 



h'^^ix) ^ M (g+^p^,^, - ^ge.pG;^^^,^}j (x, ziT'))u'^'u^'dT', (68) 
h'^^^ix) ^ M V, (^G+.pa'p' - \gapG^ (x, z{T'))u-'u^'dT'. (69) 
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In these expressions, G+ is the Lorenz gauge retarded Green's function, normahzed with a 



factor of — IGtt, following [g]. As previously mentioned, the symbol r"^ indicates that the 
range of the integral extends just short of the retarded time Tret, so that only the "tail" 
i.e., interior of the light cone) portion of the Green's function is used (see, e.g., reference 
8] for details). We define /ij^^^, rather than working with derivatives of k'^p, because /ij^^' 
is not differentiable on the worldline. (However, this non-differentiability is limited only to 
spatial derivatives of spatial components of /i^^^^, so that expressions expressions like (jl]) 
are well-defined.) A choice of retarded solution (corresponding to "no incoming radiation") 
was made in writing these equations. This choice is not necessary, and one could add an 
arbitrary smooth solution of the linearized Einstein equation to the first order in A term 
on the right side of eq. flM|) . which could then be carried through all of our calculations 
straightforwardly. However, for simplicity, we will not consider the addition of such a term. 

Our derivation of gravitational self-force to first order in A will require consideration of 
second-order metric perturbations, so we will have to carry the expansion of gab{^) somewhat 
beyond eq. (IMl) . (This should not be surprising in view of fact that our above derivation 
in section HVl of geodesic motion at zeroth order in A, required consideration of first-order 
metric perturbations.) In particular, we will need an explicit expression for the quantity 
{O'p.v)o2 appearing in the far zone expansion eq. (p^ . i.e., the term of order A^ that has the 
most singular behavior in 1/r (namely, 1/r^). 

The second-order perturbation jab satisfies the second-order Einstein equation, which 
takes the form 

G^2\j] = -G^hKh] (70) 



A=0- 



where G\fl denotes the second order Einstein tensor about the background metric g \ 
Since the 0(l/r) part of hab corresponds to the linearized Schwarzschild metric in isotropic 
coordinates (see eq. (|M|) ). it is clear that there is a particular solution to eq. (170]) of the 
form 

= — (-2W/3 + ^n^np) + O(r-i) (71) 



This can be seen from the fact that differentiation of W'^p on the worldhne 7 yields h^^^^ plus the coincidence 
limit of the integrand of (|68p. which is proportional to Raopo times the gradient of r^ct- 
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as r — ^ 0, where 

n' = x'/r (72) 



and = 0, whereas ta = Sao- (The exphcitly written term on the right side of eq. fl7T|) 
is just the O(M^) part of the Schwarzschild metric in isotropic coordinates.) The general 
solution to eq. fITOl) can then be written as 



Jab = jL + Jab , (73) 

where is a homogeneous solution of the linearized Einstein equation. We wish to compute 
the 0(l/r^) part of j^, i.e., writing 

jS = ^2^^ + 0(r-'), (74) 

we wish to compute Cab- Now, although the equations of motion to first order in A depend 
upon a choice of gauge to first order in A (see section |V]), they cannot depend upon a choice 
of gauge to second order in A, since a second order gauge transformation cannot affect the 
mass dipole moment of the background scaled metric ^^p(A = 0). (We have also verified 
by a direct, lengthy computation that second order gauge transformations do not produce 
changes in the equations of motion to first order in A.) Therefore, we are free to impose 
any (admissible) second order gauge condition on j^. It will be convenient to require that 
the Lorenz gauge condition V°(j^ — ^j^gab) = be satisfied to order The O^l/r"^) 

part of the linearized Einstein equation together with the 0(l/r^) part of the Lorenz gauge 
condition then yields 

^'^i =0 (75) 

9' (^Q,it,OA?j =0. (76) 

where Cab = Cab — -^Criab- This system of equations for -^C^u is the same system of equations 
as is satisfied by stationary solutions of the fiat spacetime linearized Einstein equation (but 
our C^^ may depend on time). The general solution of these equations is Cij = 0, Cio = 
F(t)ni + Sij{t)n^, and Cqo = 4Pi(t)n*, where Sij is antisymmetric, Sij = —Sji, where F, 
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Sij, and Pi have no spatial dependence, and where n* was defined by eq. (172|) . By a further 
second order gauge transformation (of the form = 5 fj,oF {t) / r) , we can set F{t) = 0. We 
thus obtain 

Coo(t,^,0) = 2P,(tK(^,0) (77) 
Qo{t,e,<f)) = S^,{t)n^{e,<j)) (78) 
a,{t,e,<f)) = 26^,Pk{t)n\e,<j)) , (79) 

which is of the same form as the general stationary / = 1 perturbation of Minkowski space- 
time (see, e.g., j^]), except that time dependence is allowed for Sij and Pj. As we shall 
see shortly, Sij and Pi correspond, respectively, to the spin and mass dipole moment of the 
body. 

We now may write for the metric through O(A^), 
ga/3{X; t, x') = r]af3 + B^,pj{t)x'x^ + O(r^) 



+ A ( —5^^ + hf'it, 0) + hf'it, 0)x* + MU^pit) + 0(r 



r 



JM^ 2.1 (80) 

+ A ( — {~2tatp + ^UaUp) + —Pi{t)n'5al3 + —tl^aSp)j{t)n^ 

+ ^K^pit, 9, 0) + H^p{t, 9, 0) + 0(r) j + 0{\') 

where we have introduced the unknown tensors K and H, and Sap is the antisymmetric 
tensor whose spatial components are 5"^^ and whose time components vanish, i.e., 

So^ = 0. (81) 

We now follow the strategy outlined in section |Vl We consider a smooth coordinate shift 
of the form ([52]), 

x'^ = x'^ - XAf^ix") + 0{\^), (82) 
and choose A'^ so as to make the mass dipole moment of g^^{X,to) vanish for all to- A 
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straightforward application of the coordinate transformation (182|) to the metric (IHOj) yields 



2M 



+ A (^^5a/3 + h'^'^ii 0) + /i*,%(t, 0)x' + MUaf^ii x') + 2A„,^(t, x') + 2B^,f,,{t)x'A\i, x') + 0{r^] 

/M^ 2 1 

+ ( 7^ (-2tat/3 + 3n,n^) + - [Pi(t) - MA{i, 0)] + -t(,5/3),-(t>^ 

+ ^, 0) + H^f3{i, 0, 0) + 0(r) j + 0(A3), 

(83) 

where we have "absorbed" the effects of the gauge transformation at orders A^r~^ and A^r° 
into the unknown tensors H,K. The corresponding "near zone expansion" (see eqs. fl2^ . 
(EZl), and (121) is 



2M 1 -2 /I 



+ A 



+ A^ 



+ 2A(,,^) + ii^,^ + ^ i^aSp)j,on^ + 2 [P,,o - MA,o] n^^a/?) + O j + tO (^^ 



oil]^to(l]^Po(l 



+ 0{X' 



(84) 



Notice that the indices on the left side of this equation have both a "hat" and "bar" on them 
to denote that they are components of (jab in the scaled coordinates associated with our new 
coordinates x^. By contrast, the indices on the right side have neither a "hat" nor a "bar", 
since they denote the corresponding components in the unsealed, original coordinates x'^. 
Thus, for example, Aa,i3 denotes the matrix of first partial derivatives of the x'^-components 
of Aa with respect to the x'^ coordinates^^. It also should be understood that all tensor 



Notice that the term ^(q.^^)^ arises from Taylor expanding ^(a,/3) with respect to the coordinates, so, 
in principle, the second partial derivative in this expression should be with respect to x'^ rather than x'^. 
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components appearing on the right side of eq. ( !84l) are evaluated at time to, and that Aa 
and its derivatives, as well as h'^^j^ and /ijj^^, are evaluated at = (i.e., on the worldline 7). 
Finally, the "reversals" in the roles of various terms in going from the far zone expansion of 
the metric eq. (IHOj) to the near zone expansion eq. (18^ should be noted. For example, the 
spin term ^t(^aSp)jn^ originated as a second order perturbation in the far zone, but it now 
appears as part of the background scaled metric in the near zone expansion. By contrast, 
the term BaipjX^x^ originated as part of the background metric in the far zone, but it now 
appears as a second order perturbation in the near zone expansion. 

It is easy to see from eq. fl5^ that P* — MA^ is the mass dipole moment of g^^ at time 
to- We therefore set 

A\t) = P\t)/M (85) 

for all t. Consequently, no mass dipole term will appear in our expressions below. 

Although we have "solved" for A^ in eq. f l85|) . we have not learned anything useful about 
the motion. All useful information about A^ will come from demanding that the metrics 
QabW — or, equivalently, QabW — be solutions of Einstein's equation. We may apply Ein- 
stein's equation perturbatively either via the far zone expansion or the near zone expansion. 
The resulting systems of equations are entirely equivalent, but the terms are organized very 
differently. We find it more convenient to work with the near zone expansion, and will do so 
below. We emphasize, however, that we could equally well have used the far zone perturba- 
tion expansion. We also emphasize that no new information whatsoever can be generated 
by matching the near and far zone expansions, since these expansions have already been 
fully "matched" via eqs. (El), ([27]), and (EHl). 

In the following, in order to make the notation less cumbersome, we will drop the "hat" 
on the near- zone coordinates and on the components g^^. No confusion should arise from 
this, since we will not have occassion to use the original scaled coordinates x'^ below. On the 
other hand, we will maintain the "hat" on the coordinates x^, since we will have occassion 

However, since x'^ coincides with x'^ at zeroth order in A and the A(^ci,p)'y appears at second order in A, 
we may replace the partial derivative with respect to sf^ by the partial derivative with respect to x'^ . 
However, equation ((85)) indicates clearly that solving for the displacement to center-of-mass coordinates 
is equivalent to simply determining the mass dipole moment in the original coordinates. The com- 
putations of this section may therefore be recast as simply solving enough of the second-order perturbation 
equations for the mass dipole moment — and hence the motion — to be determined. 
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to use both and below. Using this notation and setting the mass dipole terms to zero, 
eq. (IMl) becomes 



2M 1 . / 1 



+ A 



h'ap + 2A(„,^) + \k^, + ^t(.^^),n^- + 0(^^+t0 f i 



t- r •■ /1\ - /I 



+ P0 



+ 0(A3) 
(86) 



where the "dots" denote derivatives with respect to t. 

We now apply the vacuum linearized Einstein equation — up to leading order, in 1/f 
as f — s> oo — to the first order term in A appearing in eq. ( l86l) . namely 



+ ^Hc.Sp)jn' +o(^ 



t O 



(87) 



It is clear that the terms of order and t/r^ in g'^^^ cannot contribute to the linearized 



Ricci tensor to order Similarly, it is clear that the terms of order and higher 

in the background scaled metric cannot contribute to the linearized Ricci tensor to order 
so, to order we see that g^^^ satisfies the linearized Einstein equation about the 



Schwarzschild metric. It is therefore useful to expand ^^^j in tensor spherical harmonics. 

We obtain one very useful piece of information by extracting the ^ = 1, magnetic parity 
part of the linearized Ricci tensor that is even under time reversal. On account of the 
symmetries of the background Schwarzschild metric, only the ^ = 1, magnetic parity, even 
under time reversal part of the metric perturbation can contribute. Now, a general ^ = 1, 
symmetric (but not necessarily trace free) tensor field Qaj3{t,r,9,(f)) can be expanded in 
tensor spherical harmonics as (see, e.g., [30|] or [3l| equations (A16-18)) 



Qoo = Qfn' 

Q,, = Qfn^n, + Qf + Q^e,^ 
Qij = Qkn'^riiTij + n,) + Q^5^Jn'' + e\in,)n\ 



(88) 
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where the expansion coefficients Qf ,Qf ^Q'i -.Qf ,Qf ^Qf ,Qi^ ,Qf are functions of 
The three-vector index on these coefficients corresponds to the three different "m-values" 
for each C. = 1 harmonic. Thus, we see that there are a grand total of eight types of £ = 1 
tensor spherical harmonics. The six harmonics associated with labeling indices A — F are 
of electric parity, whereas the two harmonics associated with M, are of magnetic parity. 

For the metric perturbation (187|) . the "constant tensors" /i^^ and A{^a,p) are purely electric 
parity and cannot contribute. It turns out that ^KapiO, (j)) also does not contribute to the 
i = 1, magnetic parity part of the linearized Ricci tensor that is even under time reversal: 
Since Kajj is independent of i the "M" part of Kajj is odd under time reversal, whereas the 
'W part of ^Kai3{0,(j)) is pure gauge. Thus, the only term that contributes to order 
to the i = 1, magnetic parity part of of the linearized Ricci tensor that is even under time 
reversal is ^t(^aSf3)jnK Satisfaction of vacuum linearized Einstein equation at order 
requires that this term vanish. We thereby learn that 

i.e., to lowest order, the spin is parallelly propagated with respect to the background metric 
along the worldline 7. 

Having set the spin term to zero in eq. llHTI) . we may now substitute the remaining terms 
in eq. ( IHTl) into the linearized Einstein equation and set the part equal to zero. It is 
clear that we will thereby obtain relations between h^^j^, ^(a,/?), and Kais- However, these 
relations will not be of direct interest for obtaining "equations of motion" — i.e., equations 
relating A* and its time derivatives to known quantities — because the quantity of interest 
Aifi always appears in combination with the quantity A^ i, which is unrelated to the motion. 
Therefore, we shall not explicitly compute the relations arising from the linearized Einstein 
equation here. 

We now consider the information on A^ that can be obtained from the near zone second- 
order Einstein equation 

GS[^^^^] = -Gi?[^«,^«], (90) 



35 



where, from eq. (186|) . we see that 



+ o\-\+io(-A+Po' 



where we have defined 



-Da/30 = + 2^(q,/3)0 (92) 

Dafii = h^ajsi + 2A(a^/3)j + 2_BQ,j/3jy4-' . (93) 

We wish to impose the second order Einstein equation to orders and i/f^, which, as 
we shall see below, are the lowest nontrivial orders in 1/f as f ^ oo that occur. First, we 
consider G^^ij[g^^\g''^'']. The terms appearing in this quantity can be organized into terms 
of the following general forms (i) g^^^ddg^^^; (ii) dg^^^dg^^^; (iii) Tg^^^dg^^^ where T denotes 
a Christoffel symbol of the background scaled metric; (iv) TTg^^^g^^^; and (v) {dT)g^^^g^^\ 
From the form of g^^'^ together with the fact that F = 0(l/f^) and dT = 0(l/f^), it is clear 
that none of these terms can contribute to G^^,j [g^^\ g'^'^''] to order or t/f^. Therefore, 
we may treat g^"^^ as satisfying the homogeneous, vacuum linearized Einstein equation. 

We now consider the linearized Ricci tensor of the perturbation g^'^\ By inspection of 
eq. (19T1) . it might appear that terms that are 0(1) (from two partial derivatives acting on 
the "B" term) and 0(l/f) (from various terms) will arise. However, it is not difficult to 
show that the total contribution to the 0(1) and 0(l/f) terms will vanish by virtue of the 
fact that the metric gab{^ = 0) is a solution to Einstein's equation and the term proportional 
to A in eq. flHOl) satisfies the far zone linearized Einstein equation (which has already been 
imposed). It also is clear that there is no contribution of g^'^^ to the linearized Ricci tensor 
that is of order t/f^. Thus, the lowest nontrivial orders that arise in the second order 
Einstein equation are indeed 1/f^ and i/f'^, as claimed above. 

The computation of the linearized Ricci tensor to orders 1/f^ and i/f^ for the metric 
perturbation g^"^^ is quite complicated, so we will save considerable labor by focusing on 
the relevant part of the linearized Einstein equation to these orders. Our hope/expectation 
(which will be borne out by our calculation) is to obtain an equation for A\qq. Since this 
quantity is of £ = 1, electric parity type and is even under time reversal, we shall focus on 
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the i = 1, electric parity, even under time reversal part of the linearized Ricci tensor of (7^^^ 
at orders and i/f^. From eq. (l88l) . we see that the i = 1 electric parity part of the Ricci 
tensor that is 0(l/f^) and even under time reversal can be written as 

RS\i=i,+,J, = ^,RS^ (94) 
= ^ {R^n'^n.n, + i?gn,) + R^n%,) , (95) 

whereas the i = 1 part of the Ricci tensor that is 0{t/f^) and even under time reversal can 
be written as 

4^=1,+,^ = 4 W^^V., + i?f ) . (96) 



Here, in contrast to the usage of (IHHIl . Rf,Rf,Rf,R![^,Rf,R[ are "constants", i.e, they 
have no dependence on (t, f). 

We now consider the terms in g^'^^ that can contribute to these Ricci terms. The term 
BaijSjX^x^ has no £ = 1 part. Nevertheless, the i = 2 magnetic parity part of this term can, 
in effect, combine with the i = 1 magnetic parity "spin term" ^t{^aSf3)jn^ in the background 
scaled metric to produce a contribution to the linearized Ricci tensor of the correct type. 
This contribution will be proportional to 

Fi = S'^^Rkioi . (97) 

For the remaining terms in g^'^\ the "spin term" ^t(^aSi3)jn^ in the background scaled 
metric will not contribute to the relevant parts of the linearized Ricci tensor, so we may 
treat the remaining terms in g^"^^ as though they were perturbations of Schwarzschild. Thus, 
only the i = 1, electric parity, even under time reversal part of these terms can contribute. 
The remaining contributors to Rf, Rf, Rf , i?f , Rf, and R[ are 



D,oot = tDf (98) 
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and 



HooU=i,+ = Hfn' (99) 

kio\,=,^+ = kfn^m + k^ (100) 

H,,\e=i,+ = H^n'^mnj + Hf,n,) + HlnHij, (101) 

where, in eq. fl98|) . we have 

Df = Doo^ (102) 

= Aoo (103) 

Df = l{3D,\-D\,) (104) 

D[ = l{-DA + 2D\^). (105) 



(The curvature term TZajS has not appeared in the above equations because it has no £ = 1 
part.) The Df, Df, Df , Df, Hf, Hf, Hf, Hf, Kf, Kf are also "constants" in these expres- 
sions. A lengthy calculation now yields 
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Using the vacuum linearized Einstein equation R^^^j = 0, we thus obtain 6 linear equations 
for our 11 unknowns. However, in order to find "universial" behavior, we are interested in 
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relations that do not involve Hap and Kap. It can be shown that there are two such 
relations^^, namely 

- AFi - 3MDf + 2MDf - 2MDf + AMDf = , (107) 

and 

-Fi- MDf + 2MDf = . (108) 

The first equation involves and spatial derivatives of A^, and does not yield restrictions 
on the motion. However, the second equation gives desired equations of motion. Plugging 
in the definitions of the quantities appearing in eq. (11081) . we obtain 

- S''RMo^ - M{h'^t,i + '^RoM' + 2^,0.) + 2M(/i^^;'o + A^oo + Ao,^o) = 0, (109) 

where we have taken advantage of the fact (noted above) that for 00 and Oi components we 
have h^^p^ = h^ai3,'y Using the equality of mixed partials Aq^^q = v4o,oi, we obtain 

Aifio = ^S'^'-Rkioi - RojQiA^ - {h'^fl - -/ioo!i^ • (HO) 

Thus, according to the interpretation provided in section |V] above, the first order perturba- 
tive correction, Z\t), to the geodesic 7 of the background spacetime satisfies 

cPZ^ 1 ■ / ■ 1 A 

In addition, we have previously found that M and Sij are constant along 7. Taking account 
of the fact that this equation is written in Fermi normal coordinates of 7 and that = 0, 
we can rewrite this equation in a more manifestly covariant looking form as 

u'^V.iu'VbZ'') = ^/?,,/5''V-/?fe,/M^ZV-((7°"+MV)(V,/i*f -iVfc^^^^^^ . (112) 

There will be three such relations in total, because the vanishing of the mass dipole moment for all time 
implies through O(A^) in near-zone perturbation theory the vanishing of the value, time derivative, and 
second time derivative of the mass dipole at time to . The third condition should follow from the first-order 
near-zone Einstein equation, which we did not fully use. In fact, it should only be necessary to impose 
that the mass dipole have no second time derivative in order to define the motion. 
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where u"" is the tangent to 7 and u^V cSab = 0. However, it should be emphasized that this 
equation describes the perturbed motion only when the metric perturbation is in the Lorenz 
gauge (see Appendix A). 

The first term in eq. fillip (or, equivalently, in eq. (11121) ) is the "spin force" first obtained 
by Papapetrou [18]. Contributions from higher multipole moments do not appear in our 
equation because they scale to zero faster than the spin dipole moment, and thus would 
arise at higher order in A in our perturbation scheme. The second term corresponds to the 
right side of the geodesic deviation equation, and appears because the perturbed worldline is 
not (except at special points) coincident with the background worldline^^. The final term is 
the "gravitational self-force" , which is seen to take the form of a (regularized) gravitational 
force from the particle's own field. Our derivation has thus provided a rigorous justification 
of the regularization schemes that have been proposed elsewhere. 

Finally, we note that, although our analysis has many points of contact with previous 
analyses using matched asymptotic expansions, there are a number of significant differences. 
We have already noted in section HVl that our derivation of geodesic motion at zeroth order in 
A appears to differ from some other derivations j^, 1^, [l2], which do not appear to impose the 
requirement that M 7^ 0. We also have already noted that in other approaches to self-force 



, what corresponds to our scaled metric at A = is assumed to be of Schwarzschild 
form. In these other approaches, first order perturbations in the near zone expansion are 
treated as time independent, and are required to be regular on the Schwarzschild horizon. 
By contrast we make no assumptions about the time-dependence of the perturbations of 
the scaled metric beyond those that follow from our fundamental assumptions (i)-(iii) of 
section [TTl Thus, our first order perturbations are allowed to have linear dependence on 
t, and our second order perturbations can depend quadratically on t. We also make no 
assumptions about the spacetime at f < ^ and therefore impose no boundary conditions at 
small f. Finally, there is a significant difference in the manner in which the gauge conditions 
used to define the motion are imposed. In j^, O], the entire £ = 1 electric parity part of 
what corresponds to our second order near zone perturbation is set to zero without proper 



^"^ Consider a one-parameter-family wherein the initial position for a body is "moved over" smoothly with A. 
In the limit M — > 0, the body then moves on a family of geodesies of the background metric parameterized 
by A, and the perturbative description of motion should indeed be the geodesic deviation equation. 
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justification^®. By contrast, our "no mass dipole" condition applies to the background near- 
zone metric and has been justified as providing "center of mass" coordinates. 

VII. BEYOND PERTURBATION THEORY 

As already mentioned near the beginning of section |Vl the quantity in eq. (11 111) is a 
"deviation vector" defined on the background geodesic 7 that describes the first order in A 
perturbation to the motion. For any one parameter family of spacetimes Qabi^) satisfying 
the assumptions stated in section [III eq. (11 111) is therefore guaranteed to give a good approx- 
imation to the deviation from the background geodesic motion 7 as A — > 0. In other words, if 
7 is described by x*(t) = 0, then the new worldline obtained defined by = AZ' (t) is the 
correct description of motion to first order in A (when the metric perturbation is in Lorenz 
gauge) and is therefore guaranteed to be accurate at small A. However, this guarantee is of 
the form that if one wants to describe the motion accurately up to time t, then it always 
will be possible to choose A sufficiently small that a;*(t) = AZ*(t) is a good approximation 
up to time t. The guarantee is not of the form that if A is chosen to be sufficiently small, 
then a;*(t) = XZ^{t) will accurately describe the motion for all time. Indeed, for any fixed 
A > 0, it is to be expected that Z^{t) will grow large at sufficiently late times, and it is clear 
that the approximate description of motion x*(t) = AZ*(t) cannot be expected to be good 
when Z^{t) is large, since by the time the motion has deviated significantly from the original 
background geodesic 7, the motion clearly cannot be accurately described in the framework 
of being a "small correction" to 7. However, the main intended application of the first order 
corrected equations of motion is to compute motion in cases, such as inspiral, where the 
deviations from the original geodesic motion become large at late times. It is therefore clear 
that eq. (Illll) . as it stands, is useless for computing long term effects, such as inspiral. 

One possible response to the above difficulty would be to go to higher order in pertur- 
bation theory. However, it seems clear that this will not help. Although the equations of 
motion obtained from nth order perturbation theory will be more accurate than the first 
order equations, they will not have a domain of validity that is significantly larger than the 

^® Note that thepart of the £ — 1 electric parity perturbation that is relevant for obtaining equations 
of motion in 0, H is of "acceleration type" (with linear growth in r) and does not have an obvious 
interpretation in terms of a shift in the center of mass. 
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first order equations. The perturbative description at any finite order will continue to treat 
the motion as a "small deviation" from 7, and cannot be expected to describe motion accu- 
rately when the deviations are, in fact, large. In essence, by the time that the deviation from 
7 has become sufficiently large to invalidate first order perturbation theory — so that, e.g., 
the second order corrections are comparable in magnitude to the first order corrections — 
then one would expect that the (n + l)th order corrections will also be comparable to the 
nth order corrections, so nth order perturbation theory will not be accurate either. Only 
by going to all orders in perturbation theory can one expect to get an accurate, global in 
time, description of motion via perturbation theory. Of course, if one goes to all orders in 
perturbation theory, then there is little point in having done perturbation theory at all. 

Nevertheless, for a sufficiently small body of suffciently small mass, it seems clear that 
the corrections to geodesic motion should be locally small and should be locally described 
by eq. fillip . By the time these small corrections have built up and the body has deviated 
significantly from the original geodesic approximating its motion, it should then be close to 
a new geodesic, perturbing off of which should give a better approximation to the motion 
for that portion of time. One could then attempt to "patch together" such solutions to 
construct a world-line that accurately describes the motion of the particle for a longer time. 
In the limit of many such patches with small times between them, one expects the resulting 
worldline to be described by a single "self-consistent" differential equation, which should 
then well-approximate the motion as long as it remains locally close to geodesic motion. 

A simple, familiar example will help illustrate all of the above points. Consider the 
cooling of a "black body". To choose a definite problem that can be put in a framework 
similar to that considered in this paper, let us consider a body (such as a lump of hot coal) 
that is put in a box with perfect reflecting walls, but a hole of area A is cut in the this 
wall. We are interested in determining how the energy, E, of the body changes with time. 
At finite A, this is a very difficult problem, since the body will not remain in exact thermal 
equilibrium as it radiates energy through the hole. However, let us consider a one-parameter 
family of cavities where A{X) smoothly goes to zero as A — 0. When A = 0, we find that 
the energy, Eq = E{X = 0), does not change with time, and the body will remain in thermal 
equilibrium at temperature Tq for all time. When we do first order perturbation theory in 



42 



A, we will find that the first order in A correction, E^^\ to the energy satisfies^^ 

"^^^'^ -aA^'^T^ (113) 



dt 

where a is the Stefan-Boltzmann constant and A^^^ = dA/dX\x=o- Note that only the zeroth 
order temperature, Tq, enters the right side of this equation because the quantity A^^^ is 
already first order in A, so the effect of any changes in temperature would appear only to 
higher order in A. Since Tq is a constant, it is easy to integrate eq. flll3p to obtain, 

= -aA^^^T^t (114) 

Thus, first order perturbation theory approximates the behavior of E{X,t) as 

E{X, t) = Eo- XaA^^^T^t (115) 



Although this is a good approximation at early times, it is a horrible approximation at 
late times, as it predicts that the energy will go negative. If one went to second order in 
perturbation theory, one would obtain corrections to eq. f lll3p that would take into account 
the first order energy loss as well as various non-equilibrium effects. However, one would 
still be perturbing off of the non-radiating background, and the late time predictions using 
second (or any finite higher order) perturbation theory would still be very poor. 

However, there is an obvious major improvement that can be obtained by noting that if 
A is sufficiently small, then the body should remain nearly in thermal equilibrium as it loses 
energy. Therefore, although perturbation theory off of the zeroth order solution may give 
poor results at late times, first order perturbation theory off of some thermal equilibrium 
solution should give locally accurate results at all times. This suggests that if A is sufficiently 

Of course, when A becomes small compared to the typical wavelengths of the radiation (as it must as we 
let A — > 0), we enter a physical optics regime where our formulas are no longer valid. We ignore such 
effects here, just as in our above analysis of the motion of bodies in general relativity we ignored quantum 
gravity effects even though they should be important when the size of the body is smaller than the Planck 
scale. 
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small, the cooling of the body should be described by 

f - . (116) 

When supplemented with the formula, E = E{T), that relates energy to temperature when 
the body is in thermal equilibrium, this equation should provide an excellent description 
of the cooling of the body that is valid at all times. In effect, eq. flll6p takes into account 
the higher order perturbative effects (to all orders in A) associated with the cooling of the 
body, but it neglects various perturbative effects associated with the body failing to remain 
in thermal equilibrium as it cools. Equation (11161) is not an exact equation (since it does 
not take various non-equilibrium effects into account) and it is not an equation that arises 
directly from perturbation theory. Rather, it is an equation that corresponds to applying first 
order perturbation theory to a background that itself undergoes changes resulting from the 
perturbation. We will refer to such an equation as a "self-consistent perturbative equation" . 
Such equations are commonly written down for systems that can be described locally in time 
by a small deviation from a simple solution. 

How does one find a "self-consistent perturbative equation" for a given system for which 
one has derived first order perturbative equations? We do not believe that there is any 
general method for deriving a self-consistent perturbative equation. However, the following 
appear to be appropriate criteria to impose on a self-consistent perturbative equation: (1) 
It should have a well posed initial value formulation. (2) It should have the same number of 
degrees of freedom as the first order perturbative system, so that a correspondence can be 
made between initial data for the self-consistent perturbative equation and the first order 
perturbative system. (3) For corresponding initial data, the solutions to the self-consistent 
perturbative equation should be close to the corresponding solutions of the first order per- 
turbative system over the time interval for which the first order perturbative description 
should be accurate. We do not know of any reason why, for any given system, there need 
exist a self-consistent perturbative equation satisfying these criteria. In cases where a self- 
consistent perturbative equation satisfying these criteria does exist, we would not expect it 
to be unique. For example, we could modify eq. (I116P by adding suitable terms proportional 
to A"^ to the right side of this equation. 

The first order perturbative equations for the motion of a small body are that the first 
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order metric perturbation satisfies 



Wchab - 2i?U%d = -IGnMuaUt^^— , (117) 

where = corresponds to a geodesic, 7 of tlie background spacetime, and -u" is tlie tangent 
to 7. If we consider tlie retarded solution to tliis equation (wliicli automatically satisfies the 
Lorenz gauge condition), we have proven rigorously in this paper that the first order in A 
deviation of the motion from 7 satisfies 

u'^V.iu'VbZ'') = -Rtcfu'Z'^u'' - {g'^' + u''u'){V,hlf - ^VXfXu" , (US) 

with 

h'fix) = M y^^" - ^g,,G\ {x, z{r'))u^'u''dT' , (119) 

where, for simplicity, we have dropped the spin term. The MiSaTaQuWa equations 

V'Vchab - 2R\b%<i = -IQnMu^itMt) L^^liZ (120) 
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dV 



u'W = -{g""' + u^u'){Vdhlf - ^Vbhlf)u%'' , (121) 

h'fix) = M {Glba'y - \9abG\ (x, z{t')) u^'n'^'dr' , (122) 

(where one chooses the retarded solution to eq. (11201) ) are an excellent candidate for self- 
consistent perturbative equations corresponding to the above first order perturbative sys- 
tem^°. Here, M"(r) (normalized in the background metric) refers to the self-consistent motion 
z{t), rather than to a background geodesic as before. Although a proper mathematical anal- 
ysis of this integro-differential system has not been carried out, it appears plausible that our 
above criteria (l)-(3) will be satisfied by the MiSaTaQuWa equations. If so, they should 
provide a good, global in time, description of motion for problems like extreme mass ratio 
inspiral. 



The Riemann tensor term does not appear on the right side of eq. (|12ip . since in the self-consistent 
perturbative equation, the deviation from the self-consistent worldline should vanish. 
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APPENDIX A: SELF-FORCE IN AN ARBITRARY ALLOWED GAUGE 

As discussed in section |Vl the description of motion will change under first-order changes 
of gauge. Indeed, in that section, we noted that under a smooth gauge transformation, the 
description of motion changes by eq. (!56|) . However, as previously stated near the end of 
section [ITl above (see equation fl2Tl) ). the allowed coordinate freedom includes transformations 
that are not smooth at r = 0. Since such gauges may arise in practice^^, we provide here 
the expression for the first order perturbative equation of motion in an arbitrary gauge 
allowed by our assumptions. We also present the corresponding self-consistent perturbative 
equations of motion. 

As previously noted in section IVTl (see the remark below eq. fl7^ ). the equations of motion 
to first order in A depend only upon the first order gauge transformation As we have seen, 
the mass dipole moment appears at second-order in (far zone) perturbation theory, so we 
must consider the effects of first-order gauge transformations on second-order perturbations. 
This is given by g^'^^ — > g^'^^ + 5g^'^\ with 

5g^^ = {C^g^'^)at + iCy\t, (Al) 
where C denotes the Lie derivative. Equivalently, we have 

^9ab = C^cg^aH + '^^%ag^Z + ^VcVc^a) + Ve^CaV^^) + ^ c^iaV , (A2) 
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For example, the Regge- Wheeler gaugefused for perturbations of the Schwarzschild metric) is not 



smoothly related to the Lorenz gauge 19|. However, it is possible that the gauge vector is bounded 



[lot , in which case perturbations in the Regge- Wheeler gauge would satisfy our assumptions (see equa- 
tion (|A3|) ). and equations of motion could be defined. On the other hand, point particle perturbations 
expressed in radiation gauges (used for perturbations of the Kerr metric) contain a log singularity along 
a string 19], and therefore do not satisfy our assumptions. 
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where Va is the derivative operator associated with the background metric Qabi^ = 0). 
In order to satisfy the criteria on allowed gauge transformations (see equation fl2T|) ). the 
components of must be of the form 



e = F'^{t,e,<p) + 0{r) , (A3) 

i.e., C,"' cannot "blow up" at r = but it can be singular in the sense that its components 
can have direction-dependent limits. 

The mass dipole moment, P*, is one-half of the coefficient of the i = 1 part of the 
leading order, part of the second order metric perturbation, gj^ . Therefore, P* may 

be extracted from the formula, 



Svr R 



hm / gg^n^dS , (A4) 

Jr=R 



where dS is the area element on the sphere of radius R. Under the gauge transformation 
generated by we have 

SgiS = rVc^?^ + 2V'^ogi]} + rVcVo^o + V,eo V^^o + V^^oVor • (A5) 

As previously noted, for an arbitrary first-order perturbation satisfying our assumptions, we 
have 

^?So^ = -1 + ^ + 0(1), (A6) 

where M is the mass of the body. From eqs. ( lA3p . (lASp and ( 1A6[) . we see that the change 
in gl^ induced by our gauge transformation is 

4? = -^en. + X + 0(l), (A7) 
^^2d'^o9i^^+d,^od%. (A8) 

Therefore, by eq. (1A4I) . the induced change in the mass dipole moment is 

3 



5P' 

87r 



lim j (-2M^%j + Px) n'dn , (A9) 
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where dVl is the element on the unit sphere. 

Equation flA9p gives the change in the mass dipole moment induced by the possibly non- 
smooth gauge transformation generated by The corresponding change in the first order 
perturbative equation of motion is determined by the change in the smooth vector field A°- 
required to eliminate the mass dipole. Writing A"- ^ A"- + SA"-, this change is given by 

6A' = 6P'/M (AlO) 

(see eq. fl85l) ). Thus, the change ^ = + induced in the deviation vector 
describing the perturbed worldline is 



SZ' = - 

Svr r 



lim j (-2f Vij- + M'^r^x) n'dQ , (All) 

with X given by eq. (lASp . In the case where our original gauge was the Lorenz gauge, it 
follows immediately from eq. fillip that the new equation of motion for is 



'^'^^ -Ro,o'Z^ - (h^''\,o - lh^'\o''^ + S'Z' , (A12) 



where 5Z^ is given by eq. (lAlip . and where, for simplicity, we have dropped the spin term. 



We may rewrite eq. (1A12P as 



^ = -Rojo'Z^ - (h^''\,o - ^/^'"'oo'') + S'Z' + RojoSZ^ . (A13) 

Note that although eq. (1A13P provides us with the desired equation of motion in an arbitrary 
allowed gauge, the terms involving components of /i'''" must still be computed in the Lorenz 
gauge. 

Now suppose one wishes to pass to a self-consistent perturbative equation associated 
with the new choice of gauge. It is not obvious how one might wish to modify the evolution 
equations for the metric perturbations in the new gauge. (One possibility would be to 
simply use eq. fll20p and then modify the result by the addition of 2V(a^f,) but it might be 
preferable to find a new equation based on a suitable "relaxed" version of the linearized 
Einstein equation for the new gauge.) However, it appears that a natural choice of self- 
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consistent perturbative equation associated to eq. (1A13P would be 



In the case where is smooth (so that, by eq. (lAlip . we have 6Z^ = — this agrees with 



the proposal of Barack and Ori 



19|. 
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